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Abstract 

Recent results from neutrino experiments show evidence for light sterile neutrinos which do 
not have any Standard Model interactions. In this work we study the hidden interaction of sterile 
neutrinos with an ”MeV scale” gauge boson (the ngHI model) with mass Mx and leptonic coupling 
g[. By performing an analysis on the z^sHI model using the data of the MINOS neutrino experiment 
we find that the values above GxjGp = 92.4 are excluded by more than 2cr C.L., where Gp is the 
Fermi constant and Gx is the field strength of the z^sHI model. Using this model we can also probe 
other new physics scenarios. We find that the region allowed by the (5 — 2)^ discrepancy is entirely 
ruled out for Mx ^ 100 MeV. Finally, the secret interaction of sterile neutrinos has been to solve 
a conflict between the sterile neutrinos and cosmology. It is shown here that such an interaction is 
excluded by MINOS for g'^ > 1.6 x 10“^. This exclusion, however, does depend on the value of g[. 
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I. INTRODUCTION 


Most of the data collected from the neutrino oscillation experiments are in agreement 
with the 3 neutrino hypothesis [Tj. However, the observation of the reactor anomaly, which 
is a deficit of electron anti-neutrinos produced in the reactors [2], together with the results 
of the MiniBooNE experiment |3] which shows evidence for —>■ conversion, cannot 
be explained by the usual 3 neutrino scenario |1]. The most popular way to clarify these 
anomalies is to assume there exists 1 (or more) neutrino state(s) which does not have any 
weak interaction (therefore is sterile), but can mix with the active neutrinos in the Standard 
Model (SM) and change their oscillation behavior pattern. 

Although most of the anomalies seen in the neutrino sector are in favor of the sterile 
models with mass ~ 1 eV, there are conflicts between the sterile hypothesis and cosmology. 
Such light additional sterile states thermalize in the early universe through their mixing 
with the active neutrinos; therefore, we effectively have additional relativistic number of 
neutrinos which can be parametrized by AAefj. In the standard model of cosmology we 
have AAeff = 0. Massive sterile neutrinos with mass ~ 1 eV and large enough mixing angles 
to solve the reactor anomalies imply full thermalization at the early universe. This means 
that for any additional species of sterile neutrinos, we should have AAefj = 1. However, 
this is not consistent with the Big Bang Nucleosynthesis and the Planck results, which 
state ATVgff ^ O.T with 90% C.L. [^. It was recently proposed in PE] that this problem 
could be solved if the sterile neutrino state interacts with a new gauge boson X with mass 
~ a few MeV. This can easily produce a large held strength for the sterile neutrinos. In 
this way the sterile state experiences a large thermal potential which suppresses the mixing 
between the active and sterile states in the early universe. Therefore, the abundance of 
the sterile neutrinos remains small, and its impact on the Big Bang nucleosynthesis (BBN), 
Cosmic Microwave Background (CMB) and the Large Scale Structure Formation would be 
negligible, hence the sterile state can be consistent with the cosmological model. 

In this work we investigate the possibility of the sterile neutrino states interacting with a 
new gauge boson A, with mass ~ MeV, which has couplings with the sterile neutrinos and the 
charged leptons in the SM. This new interaction of the sterile neutrinos was hrst mentioned 
in [8]. The "vg Hidden Interaction^' (z^sHI) model produces a neutral current (NC) matter 
potential for the sterile states proportional to Gx, where Gx is the held strength of the new 
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interaction. The NC matter potential in the z/sHI model changes the oscillation probability 
of neutrinos and anti-neutrinos drastically. Therefore, using the data of a neutrino oscillation 
experiment such as the MINOS experiment [S], we can test the VgW. model. 

An advantage of the z/^HI model is that through it we can use the data of neutrino 
oscillation experiments to test other new physics scenarios which imply having couplings 
with a light gauge boson, such as the explanation of the {g — 2)^ discrepancy with a light 
gauge boson HDl and the secret interaction of sterile neutrinos proposed in which solves 
the tension between the sterile hypothesis and cosmology. 


II. THE FORMALISM 


We enlarge the SM with one extra species of the sterile neutrinos which do not couple 
with the SM gauge bosons, but have interactions with a new Ux{^) gauge symmetry (the 
z/gHI model). The new gauge boson couples to the sterile neutrinos and charged leptons with 
coupling constants g'^ and gl, respectively [H], where for simplicity, we have assumed equal 
coupling constants for the charged leptons. The strength of this new interaction is given by 


G 


X 


gWi 


( 1 ) 


V2 4M|.’ 

where Mx is the mass of the new gauge boson. 

The active neutrinos of the SM have charged and neutral current interactions with the 
and Z bosons. Their matter potential is therefore given by Vair) = SaeVcciz) + Vnc{'>^) = 
\/2GFNe{r){5ae “ 1/2), wliere a = e,p,r and Vcc{nc) is the charged (neutral) current 
potential of the active neutrinos. The factor Gp is the Fermi constant, while Ne{r) is 
electron number density of the earth given by the PREM model |I2|. We have assumed that 
the electron and neutron number densities are equal for our practical purposes. 

The sterile neutrinos which couple to the X boson will also have neutral current matter 
potential which is proportional to the strength held of the new interaction: 


"^sir) = -^GxNe{r) = aVNcir), 


( 2 ) 


where the dimensionless parameter a is dehned as 

Gx 


a = 


Gf 
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For a —)• 0 we recover the minimal 3 + 1 sterile neutrino model. In the minimal ”3 + 1” 
model [13] the flavor and mass eigenstates of neutrinos are related through the unitary 
(3 + 1) X (3 + 1) PMNS matrix [/://„ = The oscillation probability of neutrinos 

is described using the active-active and active-sterile mixing angles, as well as the mass 
squared differences Am^i, and Am|;^, where Am^^ = — m'j. 

The evolution of neutrinos in the VgW. model can be found by solving the following 
Schrodinger-like equation 
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where U is the 4x4 PMNS matrix |1], which is parametrized by the active-active mixing 
angles (6*12, 6*13, 6*23) as well as 3 active-sterile mixing angles (6*14, 6*24, 6*34). The matrix 

= diag^O, Amg;^, Amgj^, Am^^ j 

is the matrix of the mass squared differences. Using Eq. (|^, the matter potential matrix in 
the i^sHI model will be (after subtracting the constant VNcif) ^ I) 

= diag(Ucc(r),0,0,K(r) - UNc(r-)) 

= A/2GirAe(r)diag(^1,0,0,(5) 

The same evolution equation applies to anti-neutrinos with the replacement —)■ 

_yvs^i(j.y consider the model with a > 0. In an effective 2-neutrino scheme the 
so called MSW resonance ra happens when cos 9 = V. Since in the model the 
sterile states have nonzero matter potential, the potential would be positive in a — Vg 
system (for a > 1), which means that at energies where the resonance condition is carried 
out, converts to Ug. 

An interesting place to test the z/^HI model is the MINOS long-baseline neutrino exper¬ 
iment [H]. The MINOS experiment which has a baseline of 735 km detects both muon and 
anti-muon neutrinos, and it is one of few experiments that is both sensitive to neutrino 
and anti-neutrino oscillation probabilities. For the baseline and energy range of the MINOS 
experiment, the oscillation probabilities of the neutrinos and anti-neutrinos are very similar 
in the usual 3 neutrino scenario. However, this does not hold in the z/^HI model anymore. 
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FIG. 1. The muon neutrino and anti-neutrino survival probabilities as a function of distance over 
neutrino energy are shown in the top and bottom, respectively. The black-dashed and the red-solid 
curves correspond to the 3 and 3-1-1 neutrino models, respectively. The blue dot-dashed curves 
represent the probabilities calculated in the model for a = 150. The standard 3 neutrino 
parameters are fixed by the NUFIT best fit values m and the active-sterile mixing parameters 
are shown in the plot. The blue shaded area is the range of LjEy for the MINOS experiment [U]. 

To see how the z/^HI model affects the oscillation probability of neutrinos, we compute 
the full numerical survival probabilities for muon (anti-)neutrino in the case of the standard 
3 neutrino scenario and in the 3-1-1 and models. We show in Fig. ([^ the survival prob- 
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ability of (top) and (bottom) for the standard 3 neutrino case (black dashed curve), 
the 3 + 1 model with a = 0 (red solid curve) and the z/^HI model with a = 150 (the blue 
dot-dashed curve). To calculate the probabilities, we have hxed the 3 neutrino oscillation pa¬ 
rameters by the best £t values of NUFIT [15]: Amh = 7.5x10“®eV^, Aml^ = 2.4xlO“®eV^, 
sin^ 6 'i 2 = 0.3, sin^ 6*23 = 0.6 and sin^ 6 *i 3 = 0 .023. The values of the active-sterile mixing 
parameters in the 3+1 and models are listed in Figure 0 Bl- Comparing the 3 neu¬ 
trino case (the black dashed curve) with the 3 + 1 model (the red solid curve), we see that 
the effect of the sterile neutrino with mass squared difference = 1 eV^ is marginal 

adding only a very fast oscillation on the top of the oscillation induced by the atmospheric 
mass squared difference Amh. However, in the model we have dramatic effects both 
for neutrino and anti-neutrino survival probabilities. When the resonance condition is ful- 
hlled, we expect stronger changes for the survival probability, while for anti-neutrinos the 
changes are milder. This can be seen in the blue dot-dashed curve at the top and bottom 
of Fig. g. 


III. THE ANALYSIS 


We analyze the collected and P^ beam data in the MINOS experiment to constrain 
the a parameter in the z/^HI model. We calculate the expected number of events in each bin 
of energy by 


X {P.u,(.4v 44, AmL. AmL; a)>. . 


( 6 ) 


where 4 = sin^fly and 


is the expected number of events for no-oscillation case 
in the ith bin of energy after subtracting background [ 0 ]; while {Psur)i is the averaged 


1^,1 


—)■ survival probability in the Hh energy bin, calculated using Eq. for 


the hxed values sin^ 9u = 0.025 and sin^ 6^34 = 0 and letting the other parameters to vary. 
To analyze the full MINOS data we define the following function 




(1 + a)Af = + (1 + b)N^ - 


obs 


(^obs)2 


+ + o 1 


(7) 


where i runs over the bins of energy (23 for events and 12 for z/^ events), is the 
expected number of events dehned in Eq. g, iV® and Nf°^ are the background and observed 
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events, respectively. The 0 '°’’® = \/N°^ represents the statistical error of the observed events. 
The parameters a and b take into account the systematic uncertainties of the normalization of 
the neutrino flux and the background events respectively, with Ua = 0.016 and at = 0.2 ng. 

After combining the function for the and events and marginalizing over all 
parameters, we hnd the following best £t values: = 2.43 x 10“^ eV^, Amh = 4.35 eV^, 

^23 — 0.67, 5^4 = 0 .03, and a = 19.95, while the ratio of the value over the number of 
degrees of freedom is y^/d.o.f = 39.7/30. When we increase a from its best fit value we 
have disagreement between the z^sHI model and the MINOS data. From this we can hnd an 
upper bound for a at 2 ct C.L.: 


a < 92.4. 


( 8 ) 


Using Eq. ([^ and Eq. ([^, we can write down the coupling gl as a function of the gauge 
boson mass Mx and hxed value of a: 


9i = 


l2V2aGi 


■Mx = 5.5 X 10 


-5 


a 


M 


X 


(9) 


7 Y 92.47 ^MeV^’ 

where we have assumed the two new coupling constants in our model are related as g'^ = ■jg'i, 
in which 7 > 1. Therefore, we can use the expression above to find an exclusion region in 
the {Mx — g'l) plane. Implementing the relation above for the MINOS experiment, we arrive 
to the the black dashed curve shown in Fig. ([^. 

A light gauge boson with mass ~ MeV can be used as a unique explanation for the 3.6 a 
discrepancy between the experimental measurement and the SM prediction of the muon 
anomalous magnetic moment, {g — 2)^ [TU]. The purple shaded region in Fig. (|^ shows 
the favored 2 ct region from {g — 2)^ discrepancy. It is shown in Ref. [19] that nearly the 
entire {g — 2)^ band is excluded by various experiments if one assumes that the light gauge 
boson decays to charged leptons with branching ratio (Br) ~ 1. However, in the z^^HI model, 
the primary decay mode of the light gauge boson is into invisibles (such as the light sterile 
neutrinos) with Br ~ 1. Therefore, all the {g — 2)^ band in Fig. (|^ will be valid in the z^sHI 
model. As Fig. ([^ shows, by comparing our results on the MINOS analysis of the z^gHI 
model with the {g — 2)^ band we can exclude all the masses below Mx 100 VtT^ MeV 
with 2a C.L.. 

Another piece of information comes from neutrino trident production: the process in 
which the gAg,~ pair is produced from the scattering of off the Coulomb held of a nucleus. 
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FIG. 2. We have shown the region of interest for the model with a light gauge boson with mass 
Mx and couplings g[ and 5 ' = jg^. The result of the analysis of the model with the MINOS 
data is shown by the black dashed curve with 2a C.L. (for 7 = 30). The purple shaded region is 
the region favored by the {g — 2)^ discrepancy, while the red curve is the CCFR |T7j measurement 
of the neutrino trident cross-section m- The blue shaded region is where the tension between the 


sterile neutrino and cosmology is relieved for f{gg,Mx) = 100 and 7 = 30 (See Eq. (10) and the 
discussion after that). 


The red solid curve in Fig. ([^ represents the results of the constrains from CCFR experiment 
on measurement of the neutrino trident cross-section [IT]. As it can be seen from Fig. (|^, 
by combining the result of the CCFR experiment with our result from the MINOS analysis, 
there is only a tiny region in the {g — 2 )^ band which is allowed by all experiments. 

The sterile neutrino states with 1 eV mass have dramatic effects in cosmology due to their 
thermalization in the early universe and are disfavored by the Planck data. This tension 
could be removed if the sterile states have interactions with a light gauge boson in the so 
called secret interaction model mEo]. This will produce a temperature dependent matter 
potential for the sterile states which is Ves = [7|, where Tg is the temperature 

of the sterile sector and Ey Mx- Therefore, the oscillation of the active to sterile neutrinos 
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would be suppressed if |14ff| S> 


I 

I 2E^ 


[7]. We define the following function: 




eflf 


-fry Af ) , 


(10) 


AmlJ2E^ / 45Am^i ^Mx 

Hence, the cosmology condition in the secret interaction model would be satished if 
f{g'^,Mx) ^ 1. Similar to Eq. ([^ we can End the values of the coupling constant which 
satisfy the cosmology condition: 


g's > 


45 Amh 


M 


X 


T. 


( 11 ) 


147r2 

Assuming that the cosmology condition is satisfied for f{g'g, Mx) = 100, then using Eq. @ 
and the relation between the 2 coupling constants g'^ = the values of g'^ above 


(7: = 1.6x10-2( 


T, 


. 2 / E^/MeV . a 7 


( 12 ) 


is excluded by the MINOS analysis. Therefore, using the MINOS data we hnd that at the 
time of BBN {Ey ~ ~ 1 MeV) and for Am\^ = 1 eV^, the values of the coupling constant 

above g'g = 1.6 x 10“^ is excluded with more than 2 ct C.L. (for 7 = 30). The blue shaded 
region in Fig. (|^ shows the cosmology condition for the values mentioned above. 


IV. CONCLUSIONS 

We have investigated the possibility that the light sterile neutrinos as suggested by the 
reactor anomaly have hidden interaction with an ”MeV scale” gauge boson. In the Hidden 
Interaction (i^^HI) model, the sterile neutrinos have neutral current matter potential. There¬ 
fore, we can use the data of the neutrino experiments to constrain this model and probe 
other new physics scenarios. The held strength of this model is described by Gx- In this 
work we studied the z/^HI model using the MINOS experiment and showed that the values 
above Gx/Gp = 92.4 are excluded. 

One consequence of the z/^HI model is constraining other new physics scenarios such as 
explaining the {g — 2)^ discrepancy with a light gauge boson. We showed that using the i^sHI 
model, the {g — 2)^ region is entirely ruled out for Mx < 100a/ 7/30 MeV by the MINOS 
data. Also, the secret interaction of sterile neutrinos which is introduced in the literature 
to solve the tension between the sterile neutrinos and cosmology is excluded by MINOS for 
g'g > 1.6 X 10“^T^ for any value of Mx, where g'^ is the coupling between the sterile states and 
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the light gauge boson. We can use the data of the future neutrino oscillation experiments 
such as DUNE |2I] to further test the model and get a dehnite answer on the presence 
of the light gauge boson. 
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